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Thermodynamics provides a transparent definition of the 
free energy of density functional theory (DFT), and of its 
derivatives — the potentials, at finite temperatures T. By 
taking the T — > limit, it is shown here that both DFT and 
spin-dependent DFT (for ground states) suffer from precisely 
the same benign ambiguities: (a) charge and spin quantiza- 
tion lead to "up to a constant" indeterminacies in the poten- 
tial and the magnetic field respectively, and (b) the potential 
in empty subspaces is undetermined but irrelevant. Surpris- 
ingly, these simple facts were inaccessible within the standard 
formulation, leading to recent discussions of apparent difficul- 
ties within spin-DFT. 

PACS numbers: 71.15.Mb, 31.15.Ew, 75.20.-g 



Density functional theory (DFT) is the method of 
choice for a wide range of theoretical computations 
of electronic systems in chemistry, condensed matter 
physics, and materials science, and is also applicable to 
other many-body problems, such as classical fluids (for an 
introduction, see Ref. M). The theoretical foundations of 
DFT have repeatedly elicited discussion and reconsidera- 
tion, e.g. Refs. f^fljinost recently in the context of spin- 
dependent DFT Specifically, the first Hohenberg- 
Kohn theorem (HKI) states |8| that a given ground-state 
density distribution n(r) determines the corresponding 
external potential v(r) uniquely, up to an overall con- 
stant (an arbitrary reference energy). In contrast, in 
spin-DFT it has been shown that the densities n(r) and 
m(r) (the spin density) do not always determine the po- 
tentials v(r) and B(r) (the magnetic field), although they 
do determine the ground state \I/ uniquely, and hence 
the energy functional F[n(r), m(r)] which features in the 
second Hohenberg-Kohn theorem (HKII) is well-defined 
l&PJ- Nevertheless, many developments in spin-DFT, in- 
cluding the main practical tool — the spin-dependent 
Kohn-Sham equations |jic|| — tacitly assume a one-to- 
one correspondence between densities and potentials. 

DFT was immediately generalized to thermal ensem- 
bles by Mermin Jll]], who followed Ref. j|] closely, prov- 
ing analogues of HKI and HKII. It was soon realized 
that finite-temperature DFT does not suffer from some 
of the ambiguities of ground-state DFT (e.g., obviously 
for T > a ground-state degeneracy no longer leads 
to a one-to-many relationship between v and n), but 
only many years later it was clarified that the functional 
F[n] (or F[n, m]) can also be obtained at T > b 



Specifically, the grand potential of an electronic system 
n depends on the temperature T, the chemical poten- 
tial n, the potential v(r) and the magnetic field B(r). 
A change of representation involves a Legendre trans- 
form, e.g. replacing (i by the electron number TV as a 
free variable. In the case of the potentials v(r) and 
B(r), a functional Legendre transform must be used: 
F[n] = fl — j dr vn for (spin-independent) DFT, and 
F[n, m] = f2 — J dr vn + J dr B ■ m for spin-DFT. 

For T > 0, the simple one-to-one nature of the re- 
lationship between n(r) and v(r), or between the pair 
n, m and the pair v, B, is guaranteed by the convexity 
argument given below. For T — > 0, the ground-state for- 
malism is regained, but the degree of convexity may also 
vanish in this limit, requiring special care. The purpose 
of the present work is to establish this limit as an alter- 
native foundation for zero-temperature DFT and spin- 
DFT. The novel result is that HKI has essentially the 
same validity in ground-state spin-DFT as in DFT. 

The Hamiltonian for electrons (in a large box) is 



H=T+AW+V+B 



(1) 



with f = (-h 2 /2m)J2 a fdr->jjl l7 \7 2 'ilj ra - the kinetic en- 
ergy, W = (e 2 /2) Jdrdr' VW'VWk - r'| 
the interaction, A = 1 a parameter introduced for later 
convenience, V = J dr v(r) J2 a tpta^ra the potential 
term, and B = - J dr Y,<t,<t' B ( r ) ' V'L'w/w the mag- 
netic term |l2j [r au i is the vector of Pauli matrices, and 
in the units used the Bohr magneton /ib = !)• The term 
B is optional, and distinguishes spin-DFT from DFT. 
The free energy in the grand-canonical ensemble is 



a functional extension of standard thermodynamics 



U by 

si 



fi = -TlogS ; S = Trcxp(-(ff-/iiV)/T) , (2) 

where H is the partition function, fj, is the chemical poten- 
tial, T > is the temperature (in energy units), and the 
total particle number operator is TV = ^ CT J dr tpl^tpru . 

The partial derivatives of are the expectation value 
of the number of electrons, N — —(dfl/dfj,), and the 
entropy S = —dQ/dT. The functional derivatives give 
the density distribution n(r) = 5tt/6v(r), and for spin- 
DFT, the magnetic moment density m(r) = — <5f2/<5B(r). 

The grand potential Q is a strictly concave functional 
of v(r) and B(r). To see this, consider the operators 
X = -(H- nN)/T and Y = ~(H' - fiN)/T, where H 
and H' incorporate different potentials, denoted v and B 
for H, and v' and B' for H'. Concavity follows because 
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for any Hermitian operators X and Y and any real a 
such that X — Y ^const. and < a < 1, one has 
Tr exp(aX + (1 - a)Y) < (Tr expi") a (Tr expy) 1 "". 




FIG. 1. (a) The grand potential T) is a concave func- 
tion of the chemical potential [i (thick curve). The Legen- 
dre transform describes this curve by its family of tangents 
(dashed line): their slopes — N and their intercept with the 
vertical axis, F — SI + (J,N. (b) The resulting F(N, T) curve 
represents the Helmholtz free energy. HKII corresponds to de- 
scribing the intercept here as the minimal value of F(N) — fiN 
attainable for a given slope u, after a generalization of the 
variables \i and N to the functional variables fj, — v(r) and 
n(r). The thin lines represent the T — > limit, for which a 
range of tangents can intersect at a single point (see text). 

In thermodynamics, it is often useful to switch repre- 
sentation by employing a Legcndre transform (see Fig. 1), 
e.g. using the Helmholtz free energy F(N, T) — T) + 
fiN, where /i on the right hand side is chosen by the 
physical condition N = —(dil/dfj,), equivalent to max- 
imising the expression £1 + fiN with respect to \x. The 
DFT energy functional F[n), generalized to spin-DFT, is 
similarly introduced as a functional Legendre transform 

F[n,m] = Sl[v,B]-Jdr{n(r) (u(r)-/z) -m(r)B(r)}, (3) 

where the right hand side is to be maximized with respect 
to v and B (for DFT, m and B are simply omitted). 

The concavity of guarantees both existence and 
uniqueness of the potentials for all reasonable density dis- 
tributions, i.e. smooth distributions with n(r) > |m(r)| 
for spin-DFT, and n(r) > for DFT. The conditions of 
u-representability or A^-representability, which are nec- 
essary in the conventional approach |0,^[, do not arise 
here. Uniqueness: as the right hand side of Eq. (^) is 
concave in v and B, it cannot have more than one max- 
imum. This replaces the standard reductio ad absurdum 
argument of DFT, based on the Raylcigh-Ritz minimum 
energy principle (or a generalization thereof), and used to 
prove HKI [p|JTl| . Existence: the maximization of Eq. (j^) 
can fail only if arbitrarily large values are obtainable on 
the right hand side (e.g. for unreasonable density distri- 
butions with n < 0, very large values of v — fi can be 
taken with negligible). Any reasonable n, m set can 
be obtained as a weighted sum of distributions for which 
solutions are known to exist. As the right hand side of 
Eq. (ra) is linear in the densities, its value can never be 



larger than the weighted sum of the values of F corre- 
sponding to these known distributions — F[n, m] is con- 
vex. The relationship between potentials and densities at 
T > is thus one-to-one. It follows that the maximum 
over v and B in Eq. (m is obtained for the physical sys- 
tem for which 8Q/5v(r) — n(r) and £fi/JB(r) = m(r). 
The free-energy functional F, in the combination 

F[n,m] + J dr {n(r) (v(r) - //) - m(r) ■ B(r)} , (4) 

can be minimized with respect to the density distribu- 
tions to give back the grand potential fl. This is the 
inverse Legendre transform in thermodynamics, and cor- 
responds to HKII in DFT |^|. The derivatives here are 
SF/Sn = n — v and SF/Sm = B, as follows from Eq. (||). 
The other derivatives of F[n, m] also follow, e.g. the en- 
tropy dF/dT = d^l/dT = — S, where as usual the deriva- 
tives of F are taken at constant n(r) and m(r), while 
those of f2 are taken at constant /Lt, v(r) and B(r) (here 
F[n, m] does not depend on fi because the chemical po- 
tential appears in Q only in the combination v — fi) . 

The derivative dF/dA = dfl/dA = (W) is of spe- 
cial interest, as F\n, mj = F n i[n,m] + L dA (dF/dA), 
makes the connection [ |14| with the Kohn-Sham nonin- 
teracting system, described by F n [ which is F for A = 0. 
The expectation value of the interaction operator, (W), 
is evaluated here for that system which has the den- 
sity distributions n and m and a reduced interaction 
strength A. The major part of the interaction energy 
is given directly by the density n as the Hartree term 
E H [n\ = J drdv' e 2 n(r>(r')/2|r - r'|. Thus 

F[n, m] = F ni [n, m] + £ H [n] + F xc [n, m] (5) 

where F xc is the exchange-correlation (xc) energy, 

F xc [n, m] = J dr / xc [n, m] (r ) , (6) 
with /xc(^), the xc energy density, defined as 
/ xc [n,m](r) = f dA^ I (%A*J) > ( 7 ) 

in terms of the so-called xc-hole density 

p^,(r,r') = * (V^t,^,^) _ nAvl) t (8) 

where n a (r) = | (n(r) + am z (r)j, a = ±1 (the sum-rule 
/ dr' p™ a ,(r, r') = -8 a y holds at T = 0). The ^ and 
derivatives of Eq. (§) give 

H - v(r) = /i - w ni (r) - e<p(r) + v xc (r) 

B(r) = B ni (r)-B xc (r) (9) 

in obvious notation, retrieving the results of Kohn and 
Sham JToj ] and generalising them to spin-DFT. In prac- 
tice, one requires approximations of F xc or / xc , which 
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allow their derivatives to be taken, e.g. the local spin- 
density approximation where / xc [n, m] (r) depends only 
on n(r) and m(r). With explicit expressions for t> xc and 
B xc thus obtained (and noninteracting physics used to 
relate n,m with Uni,Bni)) Eq. (S) yields a very powerful 
and simple scheme (simpler than Hartree), for compu- 
tations of n, m and £1 of interacting electrons in given 
external fields. Eqs. (^ — 0) for F xc in terms of p xc form 
the basis for most discussions and improvements of the 
accuracy of DFT, as e.g. in the recent Ref. fLqj , which 
suggests a modified LDA for cases involving degeneracy. 

At zero temperature, and F are defined by their 
(well-behaved) T — > limits: tends to the minimal 
eigenvalue of H — fj,N, and its Legendre transform gives 
F, which at T = is the internal energy. Two issues 
require consideration: (a) the degree of concavity or con- 
vexity vanishes as T — > 0, resulting in discontinuities (in- 
determinancy) in the derivatives; and (b) interest in fully 
polarized systems, |m(r)| = n(r), becomes legitimate. 

Consider first fully polarized systems, defining the 
value of F[n, m] at |m(r)| = n(r) as the limit of its T = 
values as |m(r)| — > n(r). When B(r) and m(r) are every- 
where parallel to z, as is very often assumed in spin-DFT, 
the spin-dependent densities n^i = |(n ± m z ) and po- 
tentials v-[i = h(v =F B z ) are useful; aligning all spins in 
(say) the up direction leaves v± undetermined (for T > 0, 
this potential diverges as rz.f — > n). Indeed, the whole 
vi (r) function may vary without changing the state of the 
system, provided only that the lowest-lying spin-down 
state remains above the energies of the spin-up electrons 
present 0. However, having uj.(r) undetermined poses 
no difficulty, because that subspace is empty, ni(r) = 0. 
The analogue of this in spin-independent DFT is the sit- 
uation with n(r) = for r in a finite region (N = if this 
region covers the whole system) — the potential within 
this region is undetermined, and one may exclude it from 
the system. It is of little significance that in spin-DFT 
the density in a subspace may vanish without any diver- 
gence in the potentials (at T = 0), whereas in DFT a 
region with n(r) = may require v — > oo on (or near) its 
boundary. As shown in Ref. || , in the one-electron case 
this situation persists for non-parallel magnetic fields. 

Consider next a thermodynamic example of vanishing 
convexity (see Fig. 1): at T = 0, the relationship between 
N and /i is step-wise [jl6f — a range of /i from the ion- 
ization potential to the electron affinity corresponds to a 
single value of N, while all N between Nq and Ao + 1 
share a single value of /i. Charge is quantized, due to 
particle-number conservation, [N, H] — 0. The system 
switches between "rigid" states when /i is varied, with a 
"soft mode" occuring at each switching event — the rel- 
evant susceptibility alternates between infinite and van- 
ishing values. Rigidity might affect applications of DFT, 
e.g. indeterminancy of /i for integer A in Eq. (^J). 

In its original form, DFT deals with electron num- 
ber quantization by restricting attention to a subspace 
of fixed N. The functional derivative SF/Sn is redefined 



in terms of density variations with J dr 5n(r) — 0, and 
contains an overall arbitrary constant, corresponding to 
a choice of reference energy for the potential v(r), or to 
adding a multiple of A to H. Alternatively one may 
fix the chemical potential fj,, and define the derivatives 
through their T — > limit (barring degeneracies, this is 
the midpoint of the discontinuity) . The derivative SF/dn 
at integer J dr n = N then gives an allowable choice of 
u(r) — /x, but choices differing from it by a (small enough) 
constant give the same n(r). Similarly Stt/5v evaluated 
for a degenerate v gives the ensemble averaged n, but 
any weighted combination of the densities of the individ- 
ual ground states also corresponds to v. 

For parallel-spin DFT, the z component of the total 
spin is a second conserved, quantized quantity, [M z , H] = 
0. Again, one can limit attention to a subspace with a 
given value of M z (and of N), requiring J dr 5m z = 0. 
Shifting the magnetic field B z (r) by a constant — making 
an "arbitrary choice of the origin of reference magnetic 
field" — corresponds to adding a multiple of M z to H 
and cannot alter the state of the system as long as M z 
is fixed. The alternative of allowing M z to vary reveals 
the succession of discontinuous switching events between 
rigid spin states. Clearly, the n(r),m z (r) distributions of 
each of these states correspond to a range of potentials, 
obtainable by shifting v(r) — fi and B z (r) by two distinct 
constants (small enough not to change N or M z ). 

One may argue on physical grounds that a linear seg- 
ment in Q — a rigid state — will not occur (for |m| < n) 
without adjacent cusps. Such a linear segment in Q corre- 
sponds to a cusp or derivative-discontinuity in F. Each 
adjacent cusp in 17, representing a ground-state degen- 
eracy or soft mode, is accompanied by quantization of 
some physical quantity O which commutes with H. In 
the cases considered above, the operator N or M z can be 
constructed from the potential terms. When the operator 
O cannot be so constructed, e.g. the angular momentum 
operator L for a spherically symmetric system, a soft 
mode or cusp in Q may result, but rigidity (a cusp in F) 
can not, precisely because no combination of potentials 
is able to "push" the system in its rigid direction. 

The arguments used here are general, and apply to 
other extensions of DFT. For example, it is possible to 
construct L from the potentials of current-DFT; for den- 
sity distributions with the quantized values of electron 
number, spin and angular momentum, the potentials are 
then determined up to three terms p"7| , corresponding 
to adding multiples of each of these operators to H . A 
simpler example is a spin-singlet state with m(r) = 
in spin-DFT, where the corresponding Hamiltonian may 
conserve all three components of M. The magnetic field 
is then undetermined in both magnitude and direction: 
a small constant-B term will not change n(r) and m(r). 

Additional support for our conclusions is provided by 
the fact that generating atypical exceptions requires tun- 
ing many more parameters than are available. Indeed, all 
existing counter-examples to the spin extension of 
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HKI are either quantized or fully-polarized states. In an 
apparent exception, Ref. || identifies conditions which 
allow a magnetic field AB(r) of constant magnitude but 
space-dependent direction to be added to H without 
changing n and m (and ip), even for |m| < n. However, 
these conditions cannot be fulfilled by a many-electron 
ground state p8| . 

The thermodynamic point of view sheds new light on 
each of the five issues raised in Ref. : (i) Excited states 
are accessible to DFT and spin-DFT as thermal ensem- 
bles, for which the one-to-one nature of the potentials — 
densities relationship is guaranteed, (ii) Construction of 
accurate xc potentials from Monte Carlo or configuration 
interaction calculations can proceed in spin-DFT as in 
conventional DFT, provided one acknowledges the "up 
to a constant" nature of both the external potential and 
the magnetic field, and the irrelevance of potentials in 
empty subspaces. (hi) The presence of an excitation- 
gap need not lead to indeterminate potentials, (iv) How- 
ever, band gaps in semiconductors and in half-metallic 
ferromagnets cause a zero-temperature first-order phase 
transition: the (spin dependent) chemical potential is dis- 
continuous as a function of filling. Such derivative dis- 
continuities are of different magnitudes (and can occur 
at different density distributions) for the interacting and 
nonintcracting cases, and devising approximations to f xc 
which account for them remains a challenge. Progress 
may be achieved by identifying the bands to which the 
electron densities in Eqs. (|) belong. Obviously, other 
electronic phase transitions will pose difficulties as well, 
at least in the thermodynamic limit, (v) Finally, func- 
tional derivatives can be taken either at T — > 0+, or at 
fixed N and M, with J drdn = J drSm = 0. Chain rules, 
as in the optimized-effective potential method, then hold. 

To summarize, thermodynamic considerations followed 
by the T — ► limit can serve as an alternative foundation 
for DFT and its spin-dependent extension. At finite tem- 
peratures, the relationship between the potentials and 
the densities is guaranteed to be one-to-one by the strict 
concavity of the grand potential Q, but for T — > lin- 
ear segments may appear in J7, due to quantization of 
spin and charge. Consenquently, in the T — ground- 
state theory the densities determine the potentials (in 
non-empty subspaces) up to two spatial constants, one 
in the external potential and one in the magnetic field. 
The fact that more complicated ambiguities do not arise, 
often assumed by practitioners, has been shown here for 
the first time. One may either work at finite tempera- 
tures or fix the total electron number and the total spin 
and work within a subspace. Either option resolves all 
the difficulties of principle raised in Refs. J^Q; the re- 
maining difficulties are practical, involving the construc- 
tion of accurate xc approximations in cases of complex 
physical behavior (electronic phase transitions). 

The application of thermodynamic considerations to 
DFT may serve as an example of the unity of physics. It 
avoids some of the pitfalls encountered in the analysis of 



spin-DFT along the lines of the original Hohenberg-Kohn 
theorems. It would be interesting to compare these de- 
velopments also to a detailed analysis of spin-DFT within 
the constrained-search approach Hfl . 
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